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f^*) ■ Abstract. We show the relation between Appell's remark on the central projec- 

CNJ , tion in the dynamic of a particle and T.-Y. Thomas and Nijenhuis studies of 

• the polynomial first integrals of the geodesic flow on a space of constant curva- 

ture. Among the consequences: the space of leading terms of a quadratic first 
f""^ ■ integral of an equation q = f(q), q € M™, is isomorphic to the space of quadrilin- 

ear forms on IR n+1 having the symmetries of the Riemannian curvature tensor. 
Such a leading term also appears as a quadratic form in what we call the pro- 
jective impulsion bivector. We characterize the cases where a quadratic first 
integral comes from a Lagrangian, thus giving a geometrical interpretation and 
a converse to some recent results by Lundmark. 
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1. Introduction 



To explain what we call projective dynamics, we will explain the relation be- 
tween the dynamics of the so-called natural systems and projective dynamics. 
We compare it with the intuitive relation between the Euclidean space and the 
""P^ ■ real projective space of the same dimension. The transformation of the first 

into the second has a global aspect and a local aspect. At the global or topo- 
logical level we "add" something at infinity, namely a projective plane. This 
| gluing compactifies the space. At the local level we "forget" some structure. 

We first forget the Euclidean form to get the so-called affine space. Then we 
forget something else. We can say for example we forget the standard parallel 
transport but still keep in mind a special class of mutually indistinguishable 
parallel transports (each one corresponding to an affine chart). 

Changing the affine space into the projective space we lose the standard parame- 
trizations of the lines. In other words we lose the concept of uniform rectilinear 
motion of a point. This dynamical statement is one of the keys of projective 
dynamics. It also relates the global aspect to the local aspect: the concept 
of uniform motion does not fit with the projective line, because the point at 
infinity is never reached and the projective line is never described completely. 

A natural dynamical system is described by an equation of the form q = \7U(q), 
where q is the position of a particle in a Euclidean space, U is a function called 
the force function, and VU denotes the gradient of U. This equation can also 
model a multiparticle system. The position space, where q lives, is then the 
Cartesian product of the position spaces of each particle. 

Natural dynamics corresponds to Euclidean geometry, and we can call it Eu- 
clidean dynamics. The Euclidean form is used to define the gradient, and is a 
fundamental ingredient of the first integral of energy. If we forget the Euclidean 
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structure of the position space, we enter a larger and more basic class of systems, 
those of the form q = f(q), where f(q) is a field of forces depending only on the 
position. They are called Newton systems in [Lun], and for us they pertain to 
afhne dynamics. It is worth to mention that Jacobi [Jac] introduced what he 
called the last multiplier in the context of this class of a priori non-Hamiltonian 
systems. 

Our main claim is that there exists a projective dynamics, i.e. a dynamics where 
the position space is the real projective space, or its double covering, or more 
frequently an open subset of this double covering. We gave a description in [Alb] 
and we will give a very elementary explanation in the next section. The local 
aspects of projective dynamics will be presented in an indirect way: to forget 
some local structure, we forget the "screen" and thus consider as equivalent 
many afhne or "curved-afhne" dynamical systems. 

The global aspect of projective dynamics is the possible prolongation "after 
infinity" of the orbits in a natural system. The first example is the equation q = 
0, where the rectilinear trajectory is prolongated in the projective line. For the 
two body problem the prolongation is as expected: the hyperbolic trajectories 
continue and describe a complete ellipse of the projective space. We showed 
in [Alb] that there is also a prolongation for the motion of a particle in the 
gravitational field of any fixed repartition of mass. For the three body problem, 
the possibility of a prolongation was discussed by Chazy [Cha], who obtained 
a rather negative answer. Only some hyperbolic trajectories have a natural 
prolongation after crossing infinity. In our language the projective force field 
becomes singular at infinity. It is interesting to note that if we define the 
interaction between the bodies by an inverse cube law instead of the Newtonian 
inverse square law, this singularity disappears. There is a prolongation for all 
the unbounded orbits. 

2. Systems defined by a field of forces 

Halphen and Appell discovered the properties of the central projection in dy- 
namics (see [Alb] for the references). They considered mostly the central projec- 
tion from an affine plane in 1R 3 , where a particle moves, to another plane where 
the projected particle moves. But Appell in a remark mentioned the possibility 
of projecting on a sphere or on the two-sheeted hyperboloid in Minkowski space, 
in order to get the remarkable Kepler problem in a constant curvature space. 
We will first put AppelPs transformation at work on a larger class of systems. 

2.1. Let V be a real vector space of finite dimension n+l. A set O C 1^ \ {0} 
is said semi-conic if, for any q g ft and any A > 0, Xq g ft. We consider a 
semi-conic connected open set 51 together with 

(i) a smooth function h : Cl — > IR with the property h(Xq) = Xh(q) for any 
q g fl, X > 0, which defines the screen H = {q g fl\h(q) = 1}, 

(ii) a smooth "force field depending only on the position" / : H — > V, tangent 
to H. 

In words h is positively homogeneous of degree one. The data (i) and (ii) defines 
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a differential system that we write 

q = f + Xq. (2.1) 

The term Xq is the radial reaction. The factor A G TR is uniquely determined 
at any time in such a way that the particle remains on the screen Ti. By this 
we mean that the initial condition (q, q) being in TTi = {(q, q) G CI X V\h(q) — 
1, dh\ q (q) = 0}, the unique solution continues on TTi.. 

To see that (2.1) is a well-defined system of ordinary differential equations, we 
observe that A is determined in a unique way as a function of (q, q), as follows. 
By successive derivations of the condition h(q) = 1 we find dh\ q (q) = and 
d 2 h\ q (q,q) + dh\ q (q) = 0. The tan gency hypothesis in (ii) is dh\ q (f) — 0. By 
Euler's relation dh\ q {q) = h(q) = 1, thus X(q,q) — ~d 2 h\ q (q,q). Substituting 
this expression in (2.1) and arbitrarily extending / to CI we get an ordinary 
system of differential equations on CI x V which possesses TTi as a 2n-dimensional 
invariant submanifold. 

2.2. Change of time. We denote as well by d/dt the above dot derivative. We 
introduce another time parameter s. We denote d/ds by '. The change of time 
is defined by a — ds/dt. For any variable x we have x = dx/dt = adx/ds = ax' . 
Equation (2.1) is changed into 

a(aq')' = f\ q + Xq. (2.2) 

2.3. Central projection. We consider an equation Q" = g\q + nQ of the same 
form as (2.1), but defined on a different screen K. with equation k(Q) = 1 by a 
force field g. We project centrally Q on H, setting Q = bq. 




The change of screen is defined by the function b(q) = h(q)/k(q), positively 
homogeneous of degree on CI. We consider the motion of q G Ti defined by the 
motion of Q G K, according to Q" — <?|q + /iQ. Then q satisfies 

(bq)" = g\ bg + nbq. (2.3) 

2.4. Change of time and change of screen. Expanding and normalizing the q" 
term we compare the change of time with the change of screen: 

q" + ~q' = —f q + —q, q" + 2 T q' = -r9 b Q + {fi- -r)q- (2.4) 
a a 2 - a 1 b o b 
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If we make a = b 2 and g\b q = b~ 3 f\ q the first three terms coincide. As the 
fourth term is determined uniquely by the condition q £ TL, it is the same in 
both equations. Both equations define the same motion of q on the screen TL. 

An objection is that the field of force g\b q = b~ 3 f\ q is not tangent to the screen 
/C at the point Q = bq. Axiom (ii) is not respected. This objection is easily 
removed. The factor A in Equation (2.1) is defined in order to maintain the 
constraint h(q) = 1. If / is not tangent to the screen, we change the formula 
^(QiQ) — — d 2 h\ q (q, q) into the more general formula X(q,q) — —d 2 h\ q (q, q) — 
dh\ q (f\ q ). Now a force field / is equivalent to a force field / + "/q, with any 
7 : — > H. We call jq a radial vector field. There is a unique choice of radial 
vector field that makes the force tangent to the screen. But when we find it 
convenient we work with forces which are not tangent to the screen. 

If we extend /, previously defined on H, in 

(iii) a smooth vector field / : O — > V, positively homogeneous of degree —3, 

the solutions of Q" = /|q + on /C are sent on the solutions of q = f\ q + Xq 
on TL by central projection and change of time. 

A vector field / satisfying Axiom (iii), and defined up to the addition of a 
radial vector field, is called a projective field of force. The founding assertion of 
projective dynamics is that the data (iii) can substitute the data (i) and (ii) of 
the previous section. The screen defined by (i) appears as unessential as soon as 
we remark that another screen defines the same dynamics, except for the time 
parametrization. Let us rephrase our result. 

2.5. Definition. The restriction of the projective force field / : £1 — > V to the 
screen TL, with equation h(q) = 1, is the vector field fh -TC ^ V, tangent to TL, 
with expression fh = dh\(q A /) = / — dh(f )q. 

2.6. Definition. Given two screens TL and /C, the extended central projection is 
the application TTL — ► TIC, (q, q) ^ (Q, Q'), such that Q = bq for some b > 
and q A q = Q A Q' . 

2.7. Founding statement. A projective force field /, restricted to two screens 
TL and K. defines on each screen a tangent force field, i.e. a differential system 
of type (2.1). The extended central projection from TL to K. sends trajectory 
on trajectory, velocity vector on velocity vector, without respecting the time 
parameters. 

The example treated by Appell is the projection of the usual Kepler problem 
on the Euclidean space to the Kcpler-Serret-Killing problems on the spaces of 
constant curvature. In the first problem the screen function h is a linear form 
on V. In the second it is the square root of a non-degenerate quadratic form. 
Appell did not point out the fundamental invariant q A q = Q A Q' , that we 
called in [Alb] the projective impulsion bivector. 

3. Homogeneous form of the first integrals 

3.1. Any force field (ii) possesses a homogeneous form, the projective force field 
(iii). Any function G H defined on TTL = {(q,v) E O x V\h(q) = 1, (dh,v) = 0} 
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has a unique homogeneous form 

G(q, v) = G H (y^y, {dh, q)v - (dh, v)q) . (3.1) 

The function G is defined on Q, x V and coincides with G-u on TTt. It satisfies 

G(q,v + 1 q) = G(q,v), G(\q, X^v) = G(q,v) (3.2) 

for any 7 e IR and any A > 0. If an extended central projection maps (q, q) 
to {Q,Q'), then G(q,q) = G(Q,Q'). If G n is a first integral of System (2.1), 
if we homogenize f in a projective force field, the dynamics (2.1) obtained by 
restriction to any screen possesses G as a first integral. 

Let G(q,v) be a function satisfying (3.2). Consider the time derivative of 
G(q, q) along a solution of q = f + Xq. It defines new function G(q, v) = 
{dG/dq, v) + (dG/dv, /}. The radial vector field Xq docs not appear according 
to (dG/dv, q) = 0, which is a corollary of (3.2). The function G is indepen- 
dent of the particular screen passing through (q, q). The equation G = means 
that G is a first integral of (2.1), whatever be the screen 7i, provided the force 
field / is defined on the whole semi-conic open set £1 as a positively homoge- 
neous vector field of degree —3, in agreement with Axiom (iii). Then G satisfies 
G(q, V + fq) = G(q, v) and G(Xq, X^v) = X- 2 G{q, v). 

If the function Gn{q, v) is positively homogeneous of degree b in v, the function 
G(q, v) given by (3.1) is positively homogeneous of degree b in v and consequently 
by (3.2) positively bi-homogeneous in (q, v) with pair of degrees [b,b). The 
function G(q,v) is sum of two positively bi-homogeneous terms with respective 
pair of degrees (b — 1, b + 1) and (b — 3, b — 1). 

It is a standard deduction from such an observation that if a first integral is a 
finite sum of terms which are positively homogeneous in the velocities, one can 
decompose it as a sum of first integrals of the form Gb + Gb-2 + - • -+Gb-2k, where 
the G m are non-zero positively homogeneous functions of degree m in v, and k 
is a non-negative integer. The highest order term satisfies (dGb/dq,v) = 0, i.e. 
it is a first integral of the free motion q = Xq. 

3.2. First integrals of the free motion. Let us consider in general a first integral 
R-H of equation (2.1) where we set / = 0. The extension R constructed from 
Rrt by the process (3.1) satisfies = R(q,v) = (dR/dq,v). Consequently for 
any (g, v) in the domain of definition, any sufficiently small 7 € TR, 

R(q,v + 70) = R(q,v), R(q + iv,v) = R{q,v). (3.3) 

We assume moreover than Ru is homogeneous of degree b in the velocity v. 
As before, for any A > 0, R(Xq,v) = R(q,Xv) = X b R(q,v). If we forget the 
"sufficiently small" condition on 7 the above relations imply 

R(q, v) = R(q, V-q) = R(v, v - q) = R(v, -q), (3.4) 

a fundamental symmetry of exchange that can be seen only after homogenization 
(and that apparently was never noticed before) . A way to combine this identity 
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with the second identity (3.2) is the statement: a first integral of the free motion 
is a function of the projective impulsion q A v. We will re-state this precisely in 
the more restricted context of the polynomial first integrals. 

4. First integral which are polynomial in the velocities 

To study a system (2.1) one tries to obtain a complete base of first integrals. 
In the classical examples, it appears useless to look for first integrals which 
are not polynomial in v. It is not a priori impossible for a first integral to be 
rational in v, without being the quotient of two polynomial first integrals. But 
in most classical systems (2.1), such integral does not exist, while polynomial 
first integrals appear frequently, often in an unpredictable way. On the other 
hand, some first integrals which are real analytic in v, but not polynomial, 
appear for example in the completely repulsive three body problem. But they 
are computed through the computation of the orbits, and thus cannot give a 
priori constraints on the dynamics. 

Afhne case. We study a system q = f(q), where q moves in an open set TL of an 
affine space, i.e. we suppose that in (2.1) the screen function h is a linear form, 
and TL is the intersection of £1 C V with an affine hyperplanc. 

4.1. Proposition. Let TL C A be an open set in an affine space A. Let W be the 
tangent vector space and / : TL — > W, q *—> f(q) a locally Lipschitzian vector 
field. Let G(q, q) be a first integral of q = f(q) which is a polynomial of degree 
b in the second variable q, the coefficients being differentiable functions of q. 
Then the highest degree term Gb(q, q) is a first integral of q = 0. 

4.2. Proposition. Let TL C A be a connected open set in an affine space A. Let 
W be the vector space tangent to A. Let Rn : TL x W — ► IR, (q, q) i— ► Rnil, <?) 
be a first integral of q = 0. Assume Rn(q,v) is a polynomial in the second 
variable v <E W whose coefficients are arbitrary functions of q G TL. Then R-h is 
a polynomial in (q, v) e A x W. 

Remarks. Proposition 4.1 is standard. Proposition 4.2 improves a theorem in 
[Nij]. Nijcnhuis assumed the continuity of R-j-c, while we do not. A first integral 
is a function which is constant along any solution. This definition does not need 
any regularity. The proof of 4.2 is that, the first integral being polynomial in 
v and satisfying (3.4), it is polynomial in q, so it is a polynomial. The idea to 
use (3.4) to reprove and improve Nijcnhuis' result comes from a discussion with 
Alain Chenciner. Later we will show that the projective point of view improves 
the description of the polynomials Ru which satisfy the hypothesis of 4.2. 

Proof of 4.1. Each term of the expansion of (dG/dq,v) + (dG/dv, f) as a 
polynomial in v is zero, in particular the highest degree term (dGb/dq, v). QED 

4.3. Lemma. Let A be an m-dimensional affine space, B an n-dimcnsional affine 
space, U d A and V C B two open sets. If F : U x V — > IR, (x, y) i-> F(x,y) 
coincides onUxV with a polynomial of degree k in x whose coefficients are 
functions of y, and also coincides with a polynomial of degree I in y whose 
coefficients are functions of x, then F is the restriction to U x V of a polynomial 
in (x,y) £ Ax B. 
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Proof. Let M = (m + fc)!/m! fc! and A*" = (n + Z)!/n! 2!. Let x\ 1 . . . , xm and 
yi, . . . , yN be points in general position in U and V respectively. We mean: on 
A the unique polynomial function of degree at most k vanishing on the Xi's is 
the null polynomial. To see that there exist M points in general position in any 
open set of A, it is sufficient to see that there exists a set of points in general 
position, and that the general position property is preserved by translation and 
dilation. 

There are polynomials Pi, ... , Pm on A such that for every polynomial P of 
degree at most k, and any x e A, P(x) = J2iti P( x -i)Pi( x )- The analogous 
polynomials for B will be called Qi, ■ ■ ■ ,Qn- Then we have 

N N M 

F(x,y) =Y t F{x,y i )QM = Y t Y t F{x i ,y j )P i {x)QM- 

j=l j=l i=l 

The first equality is the hypothesis that for a given x, F(x, y) is a polynomial in 
y of degree at most I. The second is that F(x, yi) is a polynomial in x of degree 
at most k. QED 

Proof of 4.2. A first integral of q = is a function Rji such that for any fixed 
(q,v) G TL x W, the real function 7 Rh{<1 + 7V,v) of the real variable 7 is 
constant on each of its intervals of definition. Let us consider A as the affine 
hyperplane of a vector space V with equation (h,q) = 1, where h G V*. Let 
ft C V \ {0} be the semi-cone constructed on TL. We construct 

R : fl x V — ► H 

(q,v) 1 — ► R u (y-^-^,{h,q)v - {h,v)q^j. 

as in (3.1). Then R is a polynomial in v whose coefficients are functions of 
q. If (h,q) > and (h,v) > we set q = (h,q)~ 1 q, v = (h,v)~ 1 v, 7 = 
(h, q) -1 (h, v) -1 . The expression above is Rn(qo,l^ 1 { v o — qo))- Let I qv = {aq + 
(3v,a >0,f3> 0}. If I qv C ft then this expression is also i?^(fo,7 _1 (fo — <Zo)) 
because Ru is first integral of q = 0. Thus R(v, —q) = R(q, v) if I qv C fi. To 
guarantee this last hypothesis, we choose any semi-conic convex set C C ft. For 
any (q, v) G C x C, we conclude that u) = -R(w, — q). By Lemma 4.3, if C is 
open, there is a polynomial function in (x, v) which coincides with R on CxC. If 
Ci C is another open semi-conic convex set, there exists another polynomial 
coinciding with R on C\ x C\. If C fl C\ ^ {0} both polynomial coincide. So by 
connectedness i? is a polynomial function which extends to V x V, and gives by 
restriction a polynomial onAxIf coinciding with R n on W x W. QED 

Case of a general screen. The conclusions we just proved on the equation q = f 
in the affine case pass to the equation q = f + Xq on a general screen. 

4.4. Proposition. Consider System (2.1). Let G : TH — > Hi be a smooth first 
integral which at each q G TL is a polynomial on T q H of degree less than b G IN. 
In other words G = Gb + Gb-i + • • • + Go, where Gi : TH — > IR is smooth and is 
a homogeneous polynomial of degree i when restricted to any fiber T q TL. Then 
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G = Gb + Gb-2 + Gb-4 + • • • is a first integral of (2.1), and Gb is a first integral 
of q = \q. 

4.5. Proposition. Consider System (2.1) with a vector field / = 0. Let : 
TTL — > IR be a first integral. Assume R-n is at each q £ 7i a polynomial on 
T q TL of degree less than b € IN. Then Ru is the restriction to TH of polynomial 
R:VxV^JEL,(q,q)^R(q,q). 

The proof is: take / and the first integral, homogenize them, restrict to an 
affine screen, apply 4.1 and 4.2 and go back to the initial screen. The remaining 
statements in 4.4 are proved in 3.1. Here we do not care about finding the least 
possible smoothness for the screen Ji. Actually, the screen is unessential, the 
dynamics being defined, except for the time parameter, at the projective level. 

Gerard Thompson [Thn] gave a short proof of 4.2, assuming the first integral 
are sufficiently differentiable. 

5. General facts on tensors with Young tableau symmetry 

5.1. A Young tableau is given by the number r e IN of its rows, the list 
(ii, . . . , i r ) € lN r , i\ > 12 > ■ ■ ■ > i r > 1, of the lengths of its rows, and a way to 
number the boxes. For example the Young tableau (5,5,3, 1), with four rows, 
"numbered horizontally" , is represented here. 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 






14 











The shape of the same Young tableau may be given by the number c = i\ of 
columns, and the list (ji, . . . , j c ) <E IN r , ji > ji > ■ ■ ■ > j c > 1, of the lengths 
of the columns. This is convenient when we insist on the columns. The above 
tableau has 5 columns with respective lengths 4,3,3,2,2. 

But when we insist on the columns the simplest numbering is often the "vertical" 
one. So we decide that the tableau [4, 3, 3, 2, 2] is not the above one, but: 



1 


5 


8 


11 


13 


2 


6 


9 


12 


14 


3 


7 


10 






4 











The Young tableaux (5,5,3, 1) and [4,3,3,2,2] have the same shape but they 
are different, because the numbering is different. 

5.2. Given a vector space V and a Young tableau Y = (ii, . . . , i r ), we consider 
the space of multilinear forms W V* , with N = i\ + ■ ■ ■ + i r . By convention, 
we associate to the first variable the first box in the tableau, to the second the 
second box in the first row, to the i\ + 1-th the first box of the second row, etc. 
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The order corresponds to the numbering of Y. We introduce two operators S 
and A, the Young symmetrizers, acting on §Q N V* with Y. 

The operator S symmetrizes a multilinear form G @ N V* in the variables 
corresponding to boxes of the first row of Y and to the boxes of the second row 
and to boxes of the third row, etc. For example if Y = (2, 2), (S(j>)(x, y, z, t) = 
<j>(x, y, z, t) + <j>(y, x, z, t) + <j>{x, y, t, z) + <p(y, x, t, z). 

The operator A antisymmetrizes a multilinear form <fi G §Q N V* in the boxes of 
the first column of Y and in the boxes of the second column, etc. For example if 

Y = ( 2 > 2 )> (M)(x,y,z,t) = <t>(x,y,z,t)-cf)(z,y,x,t)-4>(x,t,z,y) + (f)(z,t,x,y). 

To write correctly such relations, we memorize that rows are associated to sym- 
metry, the columns to antisymmetry, and we work mentally with diagrams as 



z 


t 


X 


y 



5.3. Proposition. For any Young tableau, there is a non-zero A G IN such 
that A : <g) N V* -> W V* and S : ® W V* -» (g)* V* defined above satisfy 
AS AS = XAS and SASA = XSA. 

This result is due to [You], p. 364. A simplified proof by J. v. Neumann is 
presented in [VdW], p. 192, in [Wcl], p. 363 or in [Wc2], p. 124. These proofs 
are purely combinatorial and do not use any theoretical background. They 
concern SASA = XSA. We obtain AS AS = XAS by transposition, noticing 
that the operator *5 acts on <§Z) N V as S acts on &) N V* , i.e. by the same 
process of symmetrization, and that the same is true for A. 

5.4. Proposition. Let Y = [71, J2 ; • ■ • ,jc] De a Young tableau with c columns, 
the boxes being numbered vertically. Let N = ji + ■ ■ -+j c . The image of AS is 
the subspace of W V* generated by the elements (£1 A • • • A ^-J <g> (£1 A • • • A 
£72) ® ■ •■ ® (£1 A •■ • A£j c ), where (£1, . . . G (V*) h . A <f> G ImAS is zero if 
for any (xi,.. -,XjA G V n , 

4>(x\,X2, ■ ■ ■ ,x jl ;x 1 ,x 2 , ■ ■ -,x h ;. . .;x 1 ,x 2 , ■ ■ -,x jc ) = 0. (5.1) 



Proof. The image of S is generated by the elements £1 <g> • • • <g> ^j 1 <8> £1 ® ■ ■ ■ ® 
£j 2 <8> • • • <8> £1 <S> • • • ® £ Jc , i.e. elements having the same covector repeated along 
a same row. This proves the first claim. By the same argument (5.1) means 
that ((f>, x) — for an x G V obtained by Young symmetrization from any 
element of (g) V. This means S<p = 0. But there exists ip such that 4> — AStp. 
So ASASip = = XASip = X(f> = 0. QED 

5.5. Taking again the example Y = (2,2), numbered horizontally, the image 
of AS contains quadrilincar forms that are antisymmetric in the first and third 
boxes on one side, in the second and the fourth boxes on the other side. Fur- 
thermore, as they are obtained after antisymmetrization (1,3) of a tensor that 
was symmetric in (1, 2), the antisymmetric (1, 2, 3) part is zero. The tensor sat- 
isfies the so-called algebraic Bianchi identity. We arrived at the famous space 
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of quadrilincar forms where the Ricmann curvature tensor lives. 



1 


2 


3 


4 



The operator S establishes an isomorphism between the image of AS and the 
image of SA. In the example above, S transforms the Riemann tensor in the 
symmetrized Ricmann tensor (see e.g. [Syn], p. 54). It is a tensor that contains 
exactly the information of the Riemann tensor, but is symmetric in (1,2), in 
(3, 4), and returns zero when symmetrized in (1, 2, 3). 

The above algebraic Bianchi identity, together with the antisymmetries, char- 
acterizes the image of AS. This is a particular case of a statement true for any 
Young tableau. In Proposition 5.7 we give this general statement. In Proposition 
5.6 we give the analogous statement about the image of SA. This statements 
are useful in many applications, and are often proved in particular cases. The 
general statement appears, in a slightly weaker version than ours, at the bottom 
of p. 145 of [PeR]. This was indicated to me by Jose Maria Pozo Soler, who also 
helped me to prepare this section as well as Sections 8 and 9. As we were not 
able to find the proofs in any place, and used the results in several particular 
cases, we decided to write them. 

5.6. Proposition. Let Y = (ii,...,i r ) be a Young tableau with r rows, the 
boxes being numbered horizontally. Let N = i\ + . . . + i r . Let (si, . . . , s r ) be 
the numbers in the boxes of the first column: Sfe = i\ + ■ ■ ■ + ik-i + 1- For any 
(m, n), l<m<n<N,we denote by T™ the transposition operator acting 
on N V*, exchanging the variables m and n of an TV-linear form 4> 6 ® N V* . 
An TV-linear form </> € <§Z) N V* is element of Im SA if and only if 

(i) it is symmetric in each row: for any row k, for any (m, n), Sk < m < n < Sk+i, 

VZ4> = 4>, 

(ii) it satisfies r — 1 "algebraic Bianchi identities" : for any row fc, 1 < k < r, 

+ + t;^+ ■■■ + = o- 

5.7. Proposition. Let Y = [ji, . . . ,j c ] be a Young tableau with c columns, the 
boxes being numbered vertically. Let N = ji + . . . + j c . Let (si, . . . , s c ) be the 
numbers in the boxes of the first row: Sk = j\ + ■ ■ ■ + jk-i + 1- An N- linear 
form 4> £ ® N V* is element of Im AS if and only if 

(i) it is antisymmetric in each column: for any column k, for any (m,n), Sk < 
m<n< Sfe+i, T™(j) = -(f), 

(ii) it satisfies r — 1 "algebraic Bianchi identities" : for any column k, 1 < k < c, 

^ - r/-v - CK* CSL-i* = °- 

Remarks. The way (ii) to write the symmetrization or antisymmetrization is 
unusual but quite efficient. In each statement, the "only if part is easy, not 
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more difficult than in the case Y = (2, 2) we just treated. We now prove the 
"if part. The arguments look quite different in both proofs. 

Proof of 5.6. Consider a <j> satisfying the hypothesis and compute, for arbitrary 
(xi, . . . x r ) £ V r , 

). (5.2) 
We mean we "fill up" A<j>, in the example Y — (5, 5, 3, 1), as follows 



Xi 


X\ 


Xl 


Xl 


Xl 


X2 


X2 


Xi 


X2 


X2 


X3 


X3 


X3 






X4 











The value of A(j) is a sum with signs of values of <j>. In each term, the positions 
of the variables on each column are exchanged. We consider any of these terms, 
e.g. the one corresponding to 



^3 


Xl 


X3 


Xl 


X2 


X2 


X3 


Xl 


X2 


Xl 


X4 


X 2 


X2 






Xl 











If we apply the last of the identities (ii), this term is changed into three terms, 
none of which possess Xi in the bottom row. Such operation would work also 
if there were already some Xi in the third row. In this case several terms 
in the Bianchi identity (ii) would be equal. We would simply put them all 
in the same side of the equation, to get their common value as a function as 
the sum of the other terms, divided by an integer. We claim that repeating 
this operation, we create new terms where the Xi's are higher and higher, to 
finish with terms where they are all "at their place", in the first row. When 
necessary we put Xi in the first column using a transposition (i). We then 
place the X2's in the second row, etc. Finally we express (5.2) as a sum where 
each term is a rational number times 4>{xi, . . . , Xi; . . . ; x r , . . . , x r ). There is a 
rational number \i which does not depend neither on nor on the Xi 's, such that 
(A(f> — ji<j)){x\, . . . , Xi; . . . ; x r , . . . , x r ) = 0. This is true for any (xi, . . . , x r ) G V r 
so S{A<f> — n4>) —Q. As satisfies (i), S<f> = ii!^! . . . i r \<P- Thus is a multiple of 
SAcj) and <j) & ImSA. It remains to check that fi ^ 0. But clearly ii!^! . . . = 
X according to Proposition 5.3. QED 

5.8. Lemma. Under the hypothesis of Proposition 5.7, we have 

(iii) <f> satisfies r(r — l)/2 algebraic Bianchi identities: for any (j,k), 1 < k < 

*- Vi* - v: + i<t> — Cf*-i* = °- 

Proof. This states the "transitivity" of algebraic Bianchi identities. It is suffi- 
cient to prove the case of three columns, c = 3. As only the upper box of the 
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third column is involved, we can suppose j 3 = 1 . 



According to (ii) between Columns 2 and 3, and (i) on Column 2, we write 

(j)(xi, ...;y 1 ,...,y n ;z) = 

= -<p(x 1 , . . . ;z,y 2 , ■ ■ -,yj 2 ;yi) + <f>(xi, ■■■■,z,y 1 ,y 3 ,...,y h ;y 2 ) 

Then we expand each term in the same way using the other identity (ii) 

-<j>(x 1 ,...;z,y 2 ,...,y j2 ;y 1 ) = 

= 4>{z, x 2 ,...;x 1 ,y 2 ,...;y 1 )- 4>{z, x 1 ,x 3 ,...;x 2 ,y 2 ,...;y 1 )-\ 

(p(x-L, . . .; z,y-L,y 3 , . . . ,y j2 ;y 2 ) = 

= -<t>(z,x 2 , ...;x 1 ,y 1 ,...;y 2 ) + cj>(z, x u x 3 , . . . ; x 2 , yi, . . . ; y 2 ) 

The first terms of the right hand sides give — (f>(z, x 2 , . . . ; yi, y 2 , . . . ; x\), the 
second terms together give <f)(z, xi, . . . ; yi, y 2 , . . . ; x 2 ), etc. All together gives 
the Bianchi identity between Columns 1 and 3. QED 

5.9. Lemma. Let Y = [ji, . . . ,j c ], N = j\ + . . . + j c and G <E) N V* satisfying 
(i) and (ii) of Proposition 5.7. Let p <G V be an arbitrary vector. Then for any 
integer I, 1 < I < c, the N — Minear form <p' obtained by contracting p at the 
top of the last I columns satisfies Relations (i) and (ii) of Proposition 5.7 for the 
Young tableau Y' = [ji, . . . ,j c -i,jc-i+i -l,...,j c - 1]. 







p 


p 


p 



























Proof. Obviously (/)' satisfies Relation (i) . Relation (ii) if k < c — I is the same 
equation for <j) and </>'. Suppose k = c — I. If jk+i = 1, the column k in Y' does 
not exist and no Relation (ii) is needed. If jk+i > 2, Relation (ii) for Y can be 
written using the second top box in Column k + 1 instead of the first. It is the 
same equation for <j> and <fr' . If k > c — I, there is one term less for <f)' , but this 
term is zero, because p is repeated in the column k + 1. QED 

Proof of 5.7. Let A be the subspace of W V* defined by (i) and (ii). Note that 
(i) means A C ImA Wc know that ImAS C A and that S establishes a bijection 
between ImAS and ImSA. It is thus sufficient to prove that S : A — > ImSA is 
injective. So we want to prove that if 4> G A is such that S(f> = 0, then = 0. 
The hypothesis S4> = is equivalent to: 

#7i, gi, ... ,gj 2 ; ...; gi, ... , qj c ) = for any (q 1 , ...,q h )e V jl . 

(5.3) 
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We want to prove: (5.3) 
any n, 1 < n < ji, 



0. We first prove by induction on n that for 



X n = (f>(qi, ...,q n ;...;q 1 ,.. . , q Jk _ 1 ; q n , q 2 , ...,q jk ) = 0. 



(5.4) 



We know that X\ = 0. We write, for each i < n, the derivative of Xi with 
respect to the variable (ft, in the direction of q n . This derivative vanishes by the 
induction hypothesis Xi = 0. We get the system: 



d\X x ■ q n 
d 2 X 2 ■ q n 



= Y 1 +--- + Y 1 k ~ 1 + x n , 

= Y 2 + ■ ■ ■ + Y 2 + X n , 



(5.5) 



= 9„_iX„_i • q n = Yn-i + ■■■ + Yn-l + *n, 



where Yf = <)>{.. .;q u . . ..g^; qi, . . . ,q i - 1 ,q n ,q i+1 , . . ■ .;qi,q 2 , . . .,qj k ). 
If i > ji then Y- =0. The Bianchi identity (iii) between the columns k and I is 
X n = Y[ + ■ ■ ■ +Yj . So the sum Y{ + ■ ■ ■ + Y n _± is zero if n < ji, because it 
has only vanishing terms, and it is X n if n > j\. When we add all the equations 
(5.5), column by column, we find X n times a positive integer. Thus (5.4) holds. 

The previous argument may be used once more. We write the system d\Xi -p = 
d 2 X 2 ■ p = ■ ■ ■ = dj 1 Xj 1 - p = 0, sum all the terms and use the Bianchi identity. 
We obtain 

X = 4>{q u ... ,q h ;. .. . , q Jk _ 1 ; p, q 2 , . . .,q jk ) = 0, for any p E V. (5.6) 

We can illustrate what we just proved by the following diagram. 



9i 


9i 


9i 


9i 


9i 




92 


92 


92 


92 


93 


93 


93 






94 











= 



9i 


9i 


9i 


9i 


P 


92 


92 


92 


92 


92 


93 


93 


93 






94 











We observe now that (5.6) is the equation (5.3) written for the N — 1-linear 
form <ft' obtained contracting p e V in (ft, at the top of the last column. This 
form satisfies (i) and (ii) by Lemma 5.9. We prove our assertion (5.3) =4- <ft = 
by induction on the number of boxes N of the Young tableau. QED 



of the quadrilinear forms 

= (t>{z,t,x,y), cft(x,y,z,t) - 



5.10. Example. Consider the space B C (g) V* 
such that, for any (x,y, z,t) E V 4 , (ft{x,y,z,t) 
—<ft(y, x, z, t) and <ft{x, y, z, t) = —<j>(x, y, t, z). Suppose (ft E B and (ft(x, y, x, y) = 
for any (x,y) E V 2 . Then (ft E f^V*, i.e. (ft is completely antisymmetric. 

Proof. Consider the projector B : B — > B, (ft ift where ift is defined by 
tft{x, y, z, t) — (ft(x, y, z, t) + (ft{y, z, x, t) + <ft{z, x, y, t). As B 2 = 3B, B = kcr B © 
Im_B. We check directly that ift is antisymmetric so Im_B = f\ 4 V*. By Propo- 
sition 5.7, keri? = ImAS for the Young tableau Y = [2, 2] numbered vertically. 
We decompose (ft = (fty + ift/3 according to the decomposition of B. Finally we 
write (ft{x, y,x,y) = (fty(x,y,x,y) = 0. Thus (fty — by 5.4. QED 
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6. The Young tableau symmetry of the highest degree term 

The highest degree term of a first integral which is a polynomial in the velocities 
with coefficients depending on the position, is extended by the homogenization 
process in a polynomial (q, q) i— » R{q, q) with the fundamental properties: (i) 
for any A € R, R(M> q) = R{<li M) = ^ b R(q> 4) 7 where b £ IN is the degree, and 
(ii) for any 7 € IR, R(q, q + 79) = R(q, q), (iii) R(q + jq, q) = R(q, q). The last 
property will appear to be a consequence of the first two. 

6.1. Definition. For any b £ IN, we call V b ' b (V) the space of polynomials 
R : V x V — > IR, (q, v) 1— > R(q, v) which (i) are homogeneous of degree b in each 
variable and (ii) satisfy R(q, v + jq) = R(q, v) for any (q, v, 7) £ V x V x 1R. 

We translate these properties into properties of the "polar form" R$ of R. It is 
the 2&-linear form on V, symmetric in the first b arguments, symmetric in the 
last b arguments, such that: 

R(q, v) = R s (q, . . . , q; v, . . . , v). (6.1) 

The polar form Rg is obtained, as well-known, by a repeated differentiation of 
the polynomial R. Condition (i) fixes the number 2b of arguments. Condition 
(ii) may be written dR(q, v + r yq)/d'j = 0, which gives 

Rs(q,...,q;q,v,...,v) = 0, for any (q, v) e V x V. (6.2) 

6.2. Proposition. Let b £ IN and Rs be a 2&-linear form, symmetric in the 
first b variables, symmetric in the last b variables. The polynomial R given by 
expression (6.1) is in V b ' b (V) if and only if (6.2) is satisfied. 

The proof is what we just said. Condition (6.2) is equivalent to: the symmetriza- 
tion of Rs in its b + 1 first arguments gives zero. By Proposition 5.6, it implies 
Rs has a symmetry with Young tableau (6, b). We can state this as follows. 

6.3. Proposition. A polynomial R is in T ' (V) if and only if there exists a 
26-linear form R^ which is, for any i, 1 < i < b, antisymmetric by exchange of 
the 2i — 1-th and 2i-th variable, such that 

R(q, v) = R A (q, v; q, v; . . . ; q, v). (6.3) 

Consequently an R e V b ' b (V) satisfies R(q,v) = (-l) b R{v,q). 



q 


q 


q 


q 


q 


V 


V 


V 


V 


V 



6.4. Definition. We call V b /{V) the space of 2&-linear Rs forms satisfying 
the symmetry condition of Proposition 6.2 and the condition (6.2). We call 
Vj^ (V) the space of 2b- linear forms Ra satisfying the antisymmetry condition 
of Proposition 6.3, and the 

6.5. Algebraic Bianchi identities. For any i, 1 < i < b — 1, the antisymmetriza- 
tion of Ra in the 2i — 1-th, 2i-th and 2i + 1-th variables gives zero. 
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So, Vg (V) = ImSA, S and A being the Young symmetrizers associated to the 
Young tableau (6, b) numbered horizontally. And according to Proposition 5.7, 
Vj^ (V) = ImAS 1 , where A and 5 are the Young symmetrizers associated to the 
Young tableau with same shape but numbered vertically. This may be stated 
as follows. 

6.6. Proposition. If R G P b ' b (V), there exists a unique 

Ra e "(V) satisfying 

(6.3). The 26-linear form defines a unique ^-linear form on /\ 2 V, which is 
symmetric. It defines uniquely a polynomial Rb ■ /\ 2 V — ► IR, homogeneous 
of degree b, such that A v) — R(q,v). The polynomial Rb is completely 

determined by its values on the decomposable bivectors. 

Proof. By 5.4, R^ is linear combination of terms (£i A^^CiA^)©- ■ -<S>(£iA£2), 
with (^1,^2) € (V*) 2 . So it is linear combination of terms w ® ■■• ®w with 
u> e /\ V"* = (/\ V)*. So it is also the symmetric 6-linear form announced. 
The second part of 5.4 gives the last statement. QED 

Clearly IR © V 1A {V) © "P 2 ' 2 (V) © • • • is an algebra for the multiplication of 
polynomials. The dimension of V b ' b {V) is n(n + l) 2 (n + 2) 2 • • • (n + 6 — l) 2 (ra + 
&)/&!(&+ 1)! if dimV = n+ 1. This is a result explained in [PeR]. It was proved 
in 1954 in [FRT]. As a result about the dimension of the space of the first 
integrals of the free motion, it was stated and proved, without mention to a 
Young tableau symmetry, in [KaL] . The result in the interesting particular case 
6 = 2 was derived earlier in [Ths]. 

Through Rb the polynomial R may be interpreted as a polynomial on the Grass- 
mannian of vectorial 2-planes in V. This interpretation generalizes to p-planes, 
p > 2 (see [FuH]). More significant to us, R appears as the most general poly- 
nomial in the projective impulsion qAv. We have established that a polynomial 
first integral of the free motion, i.e. the motion with constant projective impul- 
sion, is simply a "polynomial in the projective impulsion" . 

7. Which quadratic first integrals are the Hamiltonian for a screen? 

7.1. An adapted local chart of the tangent bundle TM is obtained by differen- 
tiation from a local chart 

$:r!^IR", q^<S>(q) = ( Xl ,...,x n ), 

of the n-dimensional manifold M, where Q C M is an open set. The adapted 
local chart is 

:ra^IR 2 ", ^ ($(<?), d$|,(0) =(x u ...,x n , yi,..., y n ), 

where Tfl C TM is the inverse image of fi by the canonical projection TM — ► 
M, and q E il is the image of £ e TO by this projection. 

7.2. A second order differential equation on a manifold M is a vector field Z on 
TM such that Xi = %)i for any i, 1 < i < n, and for any adapted chart. By x we 
mean the derivative of the coordinate function Xi along the vector field Z. 
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7.3. A pre-Lagrangian for the second order differential equation Z is a function 
L : TM — > IR such that in any adapted chart (xi, . . . , i„, j/i, . . . , y n ) Lagrange 
equations 

d /dL\ dL 
dt V dyi ) dxi 

are satisfied. Here d/ dt is the derivation along the vector field Z . 

We do not require that the Lagrange equations define the system Z. We just 
require they are satisfied. The pre-Lagrangians form a vector space, which 
includes the functions Y^i=i ViVii where the j^'s are the coordinates of a closed 
1-form on M. An example of second order differential equation with many 
quadratic pre-Lagrangians is the harmonic oscillator on IR", with equations 
Xi = —Xi, 1 < i < n. The functions ±iXj — XiXj, 1 < i < j < n, are pre- 
Lagrangians. 

The following proposition is well-known in the case of Lagrangians, and the 
proof is not more complicated in the case of pre-Lagrangians. 

7.4. Proposition. Let L be a pre-Lagrangian for a second order differential 
equation Z. Given an adapted chart (x\, . . . , x n , yi, . . . , y n ), let pi = dL/dyi. 
The function YM=iPiVi does not depend on the chart and E — Y^iPiVi — L is 
a first integral of Z . 

7.5. Proposition. Let Z be a second order differential equation and G : TM — ► 
]R be a function. Using any local chart (xi, . . . , x n ) of M, we define G» : TM — > 
T*M by its local expression (xi, . . . ,x n ,yi, . . . ,y n ) ^ (xi, . . . ,x n ,pi, . . . ,p n ), 
where pi — dG/dyi. Let lo be the pre-symplectic form on TM, pull-back by G» 
of the canonical symplcctic form on T*M. Let il C M be a simply connected 
open set and TO, C TM its inverse image by the canonical projection 7r : 
TM M. The Lie derivative £zw vanishes on Tfi if and only if there exists a 
U : ft — * IR such that L = G + U o7risa pre-Lagrangian of Z on TO. 

Remark. If M is not simply connected and Cztu = on TM, then there exists 
a multivalued pre-Lagrangian for Z. It has the form G + U o n, where U is a 
multivalued function on M. 

Proof. We write 



Erasing the exact form, and using £zu = d(Z]u>) we get: w is preserved if and 
only if 




where 
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is closed. In particular, the coefficients must be independent of the y^'s. Locally 
there exists a function U (q) such that 

Ft n 

dxi 

7.6. More specific systems. We continue denoting by £ £ TM the state of the 
particle and q £ M its position, i.e. the canonical projection of £. Among the 
second order equations on M, those of the form 

D € £ = f(q), (7.1) 

where D is a torsion-free linear connection, and / is a tangent vector field on 
M, have simple properties. Here £ denotes also the velocity field along the 
trajectory. Let L be a pre-Lagrangian for this equation, which is polynomial 
in the velocity. Then both the even and the odd part of L are polynomial pre- 
Lagrangians. If L is even and of degree two, then L = T + U , where T is a 
quadratic form in the velocity £, and U only depends on the position q. There 
exists a bilinear symmetric form g on M such that 2T(£) = Lagrange's 
equations are equivalent to: Dg — and g(f, .) = dU. 

7.7. Yet more specific systems. We pass to screen dynamics. System (2.1) on 
the screen TL is of type (7.1). The connection on TL is induced by the affinc 
connection on £ and the splitting V — T q TL © [q] at any q £ TL. By [q] C V we 
mean the vector line generated by q. If D is this linear torsion-free connection, 
System (2.1) is = f(q), where £ is the velocity field along the trajectory. 
The above conclusion applies. 

7.8. Proposition. If System (2.1) on the screen Tt possesses the quadratic pre- 
Lagrangian L = T + U , the free motion q = Xq on TL possesses the quadratic 
pre-Lagrangian T . As a first integral of the free motion, T is homogenized by 
Formula (3.1) in an R £ V 2 - 2 {V). 

Which are the elements R £ V 2,2 (V) and the screens TL such that the restriction 
of R to TL is a pre-Lagrangian for the free motion? Such an R and such an TL 
are said to be compatible. We will establish some algebraic propositions in the 
next two sections before attacking this question. 



8. Maps from decomposable bivectors to decomposable bivectors 

Here we prove a result of linear algebra. The statement is more than sufficient 
to establish the result we need in Section 9. We recall that a bivector tt £ /\ 2 V 
is called decomposable if there exists (x,y) £ V 2 such that n = x Ay. We recall 
that: "7T is decomposable" <^7rA7r = 0<^>rk7r<2. 

8.1. Theorem. Let V and W be to real vector spaces with same finite dimension. 
Let Tl : A V -> A W be a linear map that sends decomposable bivectors on 
decomposable bivectors. Then either (i) there exists a non-zero <f> £ W such 
that 1Z(ir) A <f> = for any ir £ f\ 2 V or (ii) there exists a non-zero ( £ W* such 
that (ilZ(ir) = for any tt £ f\ 2 V or (iii) there exists an invertible linear map 
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B : V -> W and e = ±1 such that 7£(a; Ay) = sB{x) A B(y) for any (a;, y) G V" 2 
or (iv) dim U = dim W = 4 and there is a non-zero /i 6 /\ 4 1U and an invertible 
linear map C : V -> W* such that 1Z(x Ay) = (C(x) A C(t/))j/x. 

Remark. Conditions (iii) and (iv) arc sufficient and mutually exclusive. They 
exhaust the case of an invertible 1Z. Condition (i) is sufficient but (ii) is not and 
suggests the study of the same problem with dim W < dim V. This section is 
dedicated to the proof of 8.1. 

We set dim V = dim W = n + 1 and study the linear maps 1Z : /\ 2 V — > /\ 2 W 
that send decomposable bivectors on decomposable bivectors, i.e. 

for any (x,y) G V 2 , K(x A y) A 1Z(x A y) = 0. (8.1) 

The simplest way to obtain such an 1Z is to take a linear application B : V — > W, 
and form the unique operator B A2 : f\ 2 V —> f\ 2 W characterized by 

B A2 (x Ay) = B(x) A B(y). (8.2) 

Not all the operators satisfying (8.1) are obtained in this way. For example, if 
n = 2, i.e. dimF = 3, (8.1) is satisfied, but if rklZ = 2,1Z cannot be a B A2 . 

8.2. Lemma,. Let S : f\ 2 V x • • • x f\ 2 V — > IR, (m, . . . ,ir p ) >-> S(m, . . . ,ir p ) be 
a p-linear symmetric form. The following properties are equivalent 

(i) S(cr, . . . , a) = for any bivector a with rker < 2p — 2, 

(ii) S(iri, . . . , 7T p _2, 7T, t) = for any decomposable bivectors . . . , tt p -2, tt ; 

(iii) there exists a T e /\ 2p V* such that S(n, ...,n) = T(tt A • • • A n). 

Proof. Under condition (i) S(a, ■■■,&) = if a = X\-K\ + • • • + x p ^2^ P -2 + xn, 
where 7Ti, . . . , 7r p _ 2 and 7r are decomposable and (xi, . . . , a; p _ 2 , x) G H p_1 . As all 
the coefficients, the coefficient of x\ ■ ■ ■ x p -2% 2 in the expansion of S(a, . . . , a) 
is zero. It is a positive integer times S(iri, . . . , n p -2, 7r, it). This proves (i) 
(ii). All the coefficients in this expansion have a repeated bivector. This 
proves (ii) (i). We consider the 2p- linear form t defined by t(ui, . . . , U2 P ) = 
S(u\Au2, ■ ■ ■ , u 2p -i Au 2p ). Assuming (ii), t(u 1} u 2 , u\,u 2 , 113, 114, . . . , u 2p -2) = 0. 
By 5.10, t is antisymmetric in its first four arguments. In the same way t is 
antisymmetric in other similar sets of four arguments, and thus completely 
antisymmetric. This proves (ii) =>■ (iii). The remaining implication is standard. 

8.3. Lemma. An 1Z satisfying (8.1) defines for any p > 2 a linear application 
H(p) : /\ 2p V — > A 2P W satisfying RWfa A ■ ■ ■ A n p ) = TZ{m) A • • • A TZ(tt p ). 

Proof. Use (ii) => (iii) in Lemma 8.2, with a /\ 2p lU-valued instead of a real 
valued form. QED 

As a result, we can substitute Condition (8.1) by: 

4 

there exists a linear application 1Z^ : A V 

2 2 

for any (tt, <j) G f\ V x f\ V, TZ(n) A 1Z(a) 



f\ W such that, 
--1Z {2) (ir A a). 



(8.3) 
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This new formulation is more convenient. It gives immediately the result: 

8.4. Proposition. If TZ satisfying (8.1) is invcrtible, it maps bijectively the 
decomposable bivectors of V on the decomposable bivectors of W. 

Proof. As TZ is surjective, TZ^ reaches all the vectors of a standard basis of 
/\ 4 W, i.e. TZ {2) is surjective, and thus invcrtible. Clearly (TZ^y 1 = (^ 1 ) (2) , 
i.e. TZ~ X satisfies (8.3) and thus (8.1). QED 

8.5. The bijective case. If 1Z is bijective and n > 4, 1Z defines a bijection 
TZ-p : V(V) — > V(W), as follows. If v G V is a non-zero vector, it defines a line 
[v] C V and an n-dimensional subspace [v] A V C /\ 2 V. We claim that the 
image by 1Z of such a subspace is a [tii]AW for some non-zero w G W. This 
fact follows from the: 

8.6. Lemma. Let V be a vector space, dim V = n + 1. Let H C f\ 2 V be a, vector 
subspace such that (i) any n G H is decomposable, (ii) if a subspace K C /\ 2 
contains 7? and has the same property, then H — K. Either there exists F <zV, 
dimF = 3, such that H = /\ 2 F, or there exists a non-zero v G V such that 
i? = [u] A V. In the first case, dimi? = 3. In the second case, dimH = n. 

Proof. Let £i, ^2,-- - be a base of iJ. The two-dimensional supports of £i 
and ^ 2 have a one dimensional intersection [x] C V (if it was zero-dimensional 
£i + £ 2 would not be decomposable, if it was two-dimensional, ^ and £ 2 would 
be proportional). We write £i — x A y\ and ^2 — 1 A j/2- The support of 
£3 must intersect both supports. If it is not in [x], [£1, £2, £3] = A 2 F where 
F = [x, 2/1,2/2] and there is no £4 (the support of £4 would cut F along a line. 
A bivector r\ G [£1,^2,^3] whose support does not contain this line is such that 
£4 + f] is not decomposable). If the intersection is [x], £3 G [i]AF. For i > 3, 
£i G [x] A V by the same arguments. QED 

An invertible 1Z satisfying (8.1), if n > 4, sends an n-dimensional space of de- 
composable bivectors of V on an n-dimcnsional space of decomposable bivectors 
of W. Doing so, it defines the map TZ-p as announced. This map is bijective, the 
inverse being constructed from TZ~ X , which satisfies condition (8.1) by Proposi- 
tion 8.4. 

The bijection TZ-p sends projective lines on projective lines. To see this, observe 
that, m, u 2 , v being non-zero vectors in V, U\ Ati 2 ^ 0, the "point" [v] G V(V) 
is on the "line" [m!,m 2 ] if and only if U\ A u 2 G [v] A V. By construction of 
TZp, the line [1*1,^2] is thus sent on the line corresponding to the decomposable 
bivector 1Z(ui Au 2 ). 

By the fundamental theorem of projective geometry [Art], as dim V = dim W > 
3, there exists an invertible linear map TZp : V — > W which maps V(V) on 
V(W) as does TZp (note that if we were working with C instead of M, TZp could 
be composed with the complex conjugation). The relation between TZ and TZp 
is 

for any (u,v) G V 2 , TZ(u A v) A TZp(u) = 0. (8.4) 
Or for any (u, v) G V 2 , TZ(uAv) and TZp(u) ATZp(v) are proportional. Both term 
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define a linear map from /\ V to /\ W. The proportionality factor between 
both maps must be a constant. With the notation introducing this section, 
TZ = A(7^-p) A2 , for some non-zero AsM. 

8.7. If n = 3, the two types of maximal spaces described in Lemma 8.6 have same 
dimension. They constitute two disconnected sub- varieties of the Grassmannian 
of 3-planes in /\ 2 V. Thus if dim V — dim W = n + 1 = 4, either TZ sends the 
[u] A V's on the [v] A W's as in higher dimension, or it sends the [it] A V's on 
the /\ 2 F's (with F C W, dimF = 3.) In the first case, the conclusion is the 
same as in the case n > 4. In the second case, we compose TZ with the bijection 
* : A W -> A W*,ir h-> ttJ/j, where fj, e A ™* 

is non-zero. The composed 
map -k o TZ : f\ 2 V — > A 2 ^* i s m the first case. We have again the same 
conclusion for this map. This is the case (iv) in Theorem 8.1. 

8.8. Proposition. Under the hypothesis of Theorem 8.1, if TZ is not bijective, 
there is a non-zero decomposable bivector in (ImTZ) . 

If E is a vector space and F C E is a subspace we call F° = {£ G E* |for any ir G 
F, (£, x) =0}. This space is called the annihilator and also denoted by Ann(_F). 

Proof of 8.8 when the dimension is even. We suppose n+ 1 = dim V = dim W = 
2m. Let m G kcr7\L. We can find 7r 2 , . . . , 7r m such that -k\ A • • • A ir m ^ 0. For 
example, if rk-7ri = 2/c, we make 7T2 = • • • = iik = tti arid choose convenient 
decomposable bivectors 7Tfe+i, . . . , 7r TO . We see immediately from its definition 
that 7^ m ) = 0. Thus ImTZ is a subspace of W which does not contain any 
bivector of rank 2m. Suppose it contains a bivector ir of rank 2m — 2. There 
exists £ G 14 7 * and 77 G W* such that £ A 77 7^ and £j7r = r/j7r = 0. For 
any a G \vnTZ and any x G IR, (n + xa) A • • • A (w + xcr) = 0. In particular, 
crA7rA---A7r = 0. By contraction of £ and 77, we get {a, £ A 77) =0, using the 
hypothesis rk7r = 2m — 2. The decomposable bivector £ A 77 is in (lmTZ)° . Now 
if the maximal rank for a bivector in ImTZ is 2fc < 2m — 2, we choose n with 
this rank and ir c <E W with rank 2™ — 2 — 2k, such that 7r + ir c has rank 2to — 2. 
With the previous argument, we conclude there is a decomposable bivector in 
(ImTZ © [tt c ])° and consequently in (Imft) . QED 

To treat the odd case we need two lemmas. Recall that if n G /\ 9 W, then 
supp7r, the support of n, is determined as the intersection of the kernels of the 
£ G W* such that £j7r = 0. For any positive integer m, we set 7r Am = ir A • • • A ir 
(repeated m times). 

8.9. Lemma. If tt G f\ q W, supp?r Am C sup P 7r A P for 1 < p < m. If ir G A 2 W, 
with rank 2m, supp7r = supp7r Ap for 1 < p < m. 

Proof. £j7T Am = m(£j7r) A?r Am - 1 is proportional to (£Jtt a p) A n Am -P = p(£\Tr) A 
7t Ato_1 , so if £j7r Ap = then £j7r Am = 0. For the second statement we use 
7r = ei A e 2 H + e 2m -i A e 2m and 7r Am = e\ A • • • A e 2m . QED 

8.10. Lemma. Let e /\' /1 ^ and 4> <E W. If G supp/j, then supp(</> A 
/j) G supp/j. supp/j,, then supp(0 A /j) = supp/x © [</>]. 

Proof. The first statement is easy. For the second, consider a £ G 14 7 * with 
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£j(0A^i) = 0. Then (£,</>)// = </>A(£j/^). But there exists an r\ E W* withryj/x = 
and {<j>, rj) ^ 0. Contracting r\ in the previous relation gives (4>, ?7}£j/J = 0. 
Thus = and using again the same relation (£, (f>) = 0. So £ vanishes on 
supp/ii © [0] . QED 

Proof of 8.8 when the dimension is odd. We suppose n + 1 = dim V = dim IF = 
2m + 1. If there is a non-decomposable bivector 7Ti G ker7?,, let (eo, . . . , e n ) be 
a base of V such that tti — eo A ei + e2 A e3 + • • •. It is easy to see that the 
elements of the standard base of /\ n V a U have the form m /\tt2 A • • • A7r m , with, 
for i > 2, TTi — ejj A ei 2 , < ii < i<i < n. For example to get e\ A • • • A e„, we 
take 7r 2 = ei A e 4 , 7r 3 = e 5 A e 6 , etc. Note that this construction works because 
rk7Ti > 4. From this and the fundamental property of VS m ' stated in Lemma 
8.3, we deduce that = 0. So the bivectors in ImlZ have rank at most 

2m — 2. By the same argument as in the even case, there is a decomposable 
bivector in (Im7?.) . 

If there are only decomposable bivectors in ker 1Z we can deduce by the same ar- 
guments that rk7e< m ) < 2. ErkTZ^ = 0, we conclude with the above argument. 
If rkK^ = 1, there is a tt G IraR with tt a " 1 ^ 0, and a non-zero £ G W* such 
that £jer Am = for all a E ImK. For any such a we have A 7 r Am - 1 ) = 0. 
Then (£|ct) A 7 r Am - 1 = and £jer G supp 7 r Am - 1 = suppTr by 8.9 and 8.10 Then 
£jcr = CJtt for some ( E W* and (CJtt) A 7 r Am - 1 = 0. Thus £j7r Am = and 
(\ir = (same support). Finally £jcr = for all a E \rrrR. which gives many 
decomposable bivectors in (ImT?.) . 

The last case is rkU^ = 2. There is a tt G IraR with tt a " 1 ^ 0, a non-zero 
£ EW* with £j7r = 0, and an r\ E W* such that £ A r; ^ and £ A 77J cr A " 1 = 
for any a E \mR. Then £ A 77J (<r A 7r Am_1 ) = for any such a. Then 
r\\ ((£_|ct) A 7r Am_1 ) = 0. Then the support of the parenthesis is strictly included 
in W, so £j<7 G supp7r Am_1 by Lemma 8.10. Then £jcr = £j7r for some £ and 
the relation becomes (77 A C)j?r Am =0. If \i E /\ 2m+1 W is non-zero, this means 
(rj A C A = or £ = a£ + bn. Finally £jer = 6?7J7r and (£ A 77, cr) = for any 
cr G Imft. QED 

8.11. Proposition. If there is a non-zero decomposable bivector £Ar; G (ImT?.) C 
A 2 W* but no non-zero element ( E W* such that, for all it E f\ 2 V, (}1Z(n) = 0, 
then there is a non-zero cf> EW such that 7^(7r) A = for all 7r G A 2 V. 

Proof. We deduce from (8.1) that for all (u,v) E V 2 , 

= £ A ??J^(u A v) A A»)= (£j^(w A i>)) A (»7jft(u A v)). (8.5) 

This is the collinearity of two elements of IF. To see how this collinearity may 
occur, we take the value of the above expression on (a, (3) E (IF*) 2 . We get 

(TZ(u A v), £ A a)(R{u A v), 77 A 0) = (H(u A v), £ A /3){H(u Av),r]Aa). 

The relation says that a polynomial in the coordinates of the vectors u, v, 
a, j3 possesses two factorizations. As the degrees are the same, if the non-zero 
polynomial (1Z(uAv),1;A/3) divides the polynomial (1Z(uAv), T]A(3), the quotient 
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is a real number A. Then (rj — A£)j7£(it A v) = for all (u, v) G V 2 , which is 
excluded by the hypothesis. 

But (K(«Ad),^ A/3) cannot divide the other factor. So it cannot be irreducible 
as a polynomial in the coordinates of u, v, (3. We claim that the factorization is 

2 

(1Z(u A f), £ A /3) = {u>, u A v)((p, (3), with ue/\V*, cj>&W. (8.6) 

The proof is as follows. The total degree in the coordinates of (3 is one for 
the product, so it should be one for a factor, zero for the other factor. As the 
product is homogeneous, the factors are homogeneous. The same is true for 
the coordinates of u and of v. So the factorization is: linear times bilinear. 
But u and v cannot be separated, because 1Z is antisymmetric in (u,v). So 
there exist u> and 4> as claimed. We have £,]7l(u A v) = (uj,u A v)(j>. We have 
assumed that this quantity is not identically zero. Putting it in (8.5) we find 
(r]}K{u A v)) A 4> = 0. From (8.6) we deduce that any (3 G ker (f> is such that 
£ A (3 G (Im7?.) . So we have Relation (8.5) with (3 instead of ij, and we deduce 
(PiTl(u A v)) A 4> = 0. Finally TZ(u A v) A <f> = 0. QED 

9. From decomposable bivectors to decomposable 2-forms 

Theorem 8.1 gives enough results on an 

ft : A V -> A W satisfying (8.1) 
to deduce more complete results on the particular case where W = V* , and 
where TZ is associated to a quadrilinear form in V^ 2 V*, i.e. a quadrilinear 
form with the symmetries of the Ricmannian curvature tensor (see 6.4 and 
5.5). Let : V A — » 1R, (u, v,w, x) i— > J?^t(u, a;) be such a form. It 
satisfies Ra{u, v; w, x) = — Ra{v , u\ w , x) = —R^(u,v;x,w) and the algebraic 
Bianchi identity Ra( u , v; w, x) + Ra(v, w; u, x) + Ra{w, u; v, x) = 0. It defines 
a linear map 

ft : A V -» A V*. This map is symmetric: (lZ(u A v ), w A x) = 
(1Z(w Ax),uAv). We write 1Z(u A v) = Ra( u , v ', ■, •)■ We define: keri?^ = {u e 
V, R A (u, .; ., .) = 0}. We write Condition (8.1) 

Ra(u, v; ., .) A Ra{u,v; ., .) = for any (u, v) eV 2 . (9.1) 

9.1. Proposition. Let V be a real vector space, G T 5 ^' 2 !^*. If -R^i satisfies 
(9.1), then the induced quadrilinear form Ra ■ (V/ ker Ra) 4 —> IR satisfies (9.1). 

9.2. Theorem. Let V be a real vector space, n + 1 = dim V > 3. Let Ra G 

and ^ : A 2 ^ ^ A 2 V* associated to it. Suppose R A satisfies (9.1) 
and keri?_4 = {0}. Then either (i) there exist an invertible symmetric linear 
map B : V — > V* and e = ±1 such that 7?. = eB A2 , or (ii) there exist a non- 
zero G V* and a non-degenerate quadratic form g defined on ker <f) such that 
Ra( u i v> u > v ) = s( ( / , J ( u A w)) . 

9.3. Remark. The second member of the last formula is well-defined: 4>\{u A 
v) = <p(u)v — 4>{v)u G ker0. The formula uniquely defines Ra according to 
Proposition 5.4. We will need the complete formula for Ra- Let g be the 
bilinear form associated to g. The complete formula is Ra{u 1 v; w, x) — g{cj>\ (uA 
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v), cf)\ (wAx)) , simply because this formula satisfies the algebraic Bianchi identity, 
as easily seen using an arbitrary extension of g to V. 

Proof of 9.1. The support of Ra, seen as a tensor of ® V*, is (kcri?^) . 
The support of Ra(u, v; ., .) G /\ 2 V* is included in (keri?^) . If Ra(u, v; ., .) 
is decomposable as an element of /\ 2 V* it is decomposable as an element of 
A 2 (keri?^)°. QED 

Proof of 9.2. We begin with the case where 1Z is invertible and n > 4. By 
Theorem 8.1, there exists a linear map B : V — > V* and an s = ±1 such that 
7Z = eB A2 . We choose e = 1. The other choice of sign is similar. Let 6 be 
the bilinear form associated to B, i.e. b(u,.) = B{u). Then Ra(u,v;w,x) = 
b(u, w)b(v, x) — b(u, x)b(v, w). We write 

Ra{u, v; w, x)+Ra{v, w; u, x)+Ra{w, u; v, x) = b(u, x)(b(w, v)—b(v, w))+pcrm. 

By the Bianchi identity this expression vanishes. It vanishes for any x G V so 

(b(w, v) — b(v, w))B(u) + (b(u, w) — b(w, u))B(v) + (b(v, u) -b(u,v))B(w) = 0. 

We claim that b(w,v) = b(v,w) for any (v,w) with v A w ^ 0, and thus b is 
symmetric. We take u such that uAvAw ^ 0. By invertibility of B, B(u), B{v) 
and B(w) are independent. The identity appears as a zero linear combination 
of these three covectors. So the coefficients are zero. 

We continue with an invertible 1Z, considering now n = 3. By 8.7, the situation 
is the same as if n > 4, except for the case where 1Z is the composition of an 
TV : f\ 2 V — > /\ 2 V with the isomorphism * : f\ 2 V — > f\ 2 V* . In this case there 
is a C : V — > V such that 1Z' = ±C A2 . We claim that such a composition does 
not give an R^ which satisfies the Bianchi identity. Identifying /\ 4 V with IR 
using *, we write Ra( u , v,w,x) = u A v A C(w) A C{x). The Bianchi identity is 

mAdA C{w) + u A C(v) Aw + C(u) A v A w = 0. 

Writing this identity on base vectors e ,ei,e2,e 3 we deduce easily that the 
matrix of C is zero. So this case is excluded, and the proof is complete in the 
case of an invertible 1Z. 

We suppose TZ is not invertible. As ker R^ = {0} excludes Case (ii) of Theorem 
8.1, we are in Case (i). There is a non-zero <fi G V* such that Ra{u 1 v; ., .) A(f> = 
for all (u, v) G V 2 . It implies Ra(u,v;w,x) — for any (w,x) G (ker</>) 2 . 

An Ra G V^iV) associates to a u G V the quadratic form v ^ Ra(u,v;u,v). 
In the present case, we have for any w G ker</>, Ra(u + w,v;u + w,v) — 
Ra{u,v;u,v), provided v G kcr0. Then any u G V such that (</>, u) = 1 defines 
the same quadratic form g on ker0. If {(f), u) ^ 1, the quadratic form defined 
by u on ker^> is still proportional to g. If g was degenerate, we would have a 
non-zero v G ker</> such that Ra(u, v;u, .) = for any u. By Lemma 5.4 this 
would imply Ra(-,v; ., .) = 0, which would contradict keri?^ = {0}. The last 
statement to prove is the formula for Ra- We have 

g{<f)\ (uAv)) = Ra( w , 4>\ ("Ad); w, 4>\(u Av)) where w G V satisfies (</>, w) = 1. 
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Suppose (4>,u) ^ and set w = u/((f>,u). The second member is R A (u, v; u, v) 
as claimed. If {4>,u) = we rather set w = v/{(f>,v) to get the same result. If 
{<jy, u) = (0, v) = the equality in 9.2 is = 0. QED 

The description in Theorem 9.2 is not completely satisfying because we do not 
see if one can pass continuously from Case (i) to Case (ii). For information, we 
give without proof the following result. 

9.4. Proposition. Let V be a real vector space, dimF = n+ 1 > 3. We choose a 
non-zero element of /\ n+1 V* and define the isomorphism 7k- : /\ 2 V — > /\™ 1 V* 
associated to it. Let G : V* — > V, satisfying t G = G. We construct an 
ft : A 2 V -» A 2 V* conjugating by * the map G An_1 : A 11-1 V* -> /\ n ^ V. We 
set R A {u,v; w,x) = (1Z(u A v), w A or). Then i?^ G P^' (V) and satisfies (9.1). 
We have keri?^ = {0} if and only if rkG > n. If G has rank n + 1, we are in 
Case (i) of Theorem 9.2, with J3 = G _1 . If G has rank n, we are in Case (ii), 
with [4>] = ker G. 

10. Free motion with quadratic pre-Lagrangian 

Using the results of Section 9 we can solve the problem raised at the end of 
Section 7. We first derive the equation for a compatible pair (R,7i). The bi- 
quadratic polynomial R : V 2 — > H, (q, v) _R(q, v) defines a quadratic form 
on the tangent space T q H at any point q e H. The compatibility condition is 
simply that this form is invariant by parallel transport. The parallel transport 
of a vector w € T q TL along the curve with tangent velocity q s T q TL is defined 
by the equation w = Xq, where A is a real function. We use the antisymmetric 
polarization R A of i?, which is such that R(q, w) = Ra(q, w; q, w). We write 

We arrive at a first form of the compatibility condition: 

for any q e H, u G T q H, w G T q H, R A {q, w; u, w) = 0. (10.1) 

According to 5.9 we can apply 5.4 to the case of the trilinear form R A {q, ■', ■, •)■ 
Thus (10.1) is equivalent to 

for any q G H, u G T q H, v G T q H, w G T q H, R A {q, v; u, w) = 0. (10.2) 

The last expression depends on q and v through the bivector q A v. We can 
change the hypothesis v G T q TL into v G V. Thinking of R A (q,v; ., .) as a 2- 
form in the missing arguments, and dh as the 1-form, differential of the screen 
function h, the compatibility condition (10.2) reads 

for any q G H,v G V, dh\ q A R A (q, v; ., .) = 0. (10.3) 

In particular, R A satisfies (9.1). The two following results correspond to Propo- 
sition 9.1. 



24 



10.1. Lemma. Let (Ra,TL) be a compatible pair. Let h(q) = 1 be the equation 
of the screen TL. If k E keri?^, then (dh\ q , k) = for any g E TL. 

Proof. If q ^ keri?_4, there exists acel/ such that Ra(<1, v \ ■, •) 7^ 0. If g G W, 
we choose (10.3) as the compatibility equation and contract k E keri?^. It 
gives (dh\ q ,k) = 0. Suppose q E TL f~l kcri?_4. By continuity we should have 
(dh\ q ,k) = 0, and in particular (c?/i| g ,g) = /i(g) = 0. This is impossible, so 
nr\kcrR A = 0. QED. 

10.2. Proposition. Let R A G 7^' 2 (V) not equal to zero, fl^ G V A a {V/ kcr i?^) 
be the quadrilinear form induced on the quotient space V/kerR^. Then any 
screen TL° of V/keri?^ compatible with R A gives by pull-back a "cylindric" 
screen TL compatible with Ra- Reciprocally every screen compatible with R A 
is contained in such a cylindric screen. 

Proof. We have T q oTL° = T q TL/ ker Ra for any q projecting on q° by the canon- 
ical projection V — > V/keri?,A. Denote by w° and w° the respective canonical 
projections of v and u>. We have R A (q°, v°; w°, w°) — RA.(q,v;w,w). If the com- 
patibility condition (10.1) is satisfied for (R A ,H°), it is satisfied for (Ra,H), 
and reciprocally. QED 

10.3. Proposition. Suppose (Ra,7~1) is a compatible pair with kcri?.4 = {0}. 
Then either (i) there exists a non-degenerate quadratic form g on V and a 
A G 1R such that the screen H is included in the quadric of equation g(q) = 1 
and Ra{q, v i <1, v ) = ^g( v ) for any q eH and v G T q H, or (ii) there exists a non- 
zero <f) G V* , a, non-degenerate quadratic form g defined on ker <j> and a A G 1R 
such that the screen TL is included in the hyperplane of equation {(f>, q) = 1 and 
RA{q, v :Q, v ) = Ag(u) for any q e H and w G T q H = ker 0, or (iii) n = 1 i.e. 
dimV" = 2. 

Proof. As (10.3) implies (9.1), we are in case (i) or (ii) of Theorem 9.2. In case 
(i) Equation (10.3) becomes dh\ q A B(q) A B(v) = 0. This is true for any v so if 
dimV > 3 we have dh\ q AB(q) = for any q E H. We set 5(g) = (B(q),q). The 
equation is c?/i| g A c?<?| 9 = 0. The screen being a level hypersurface of the screen 
function h, it coincides according to the equation with a level hypersurface 
of g. We compute Ra{q, v ', v ) = (R-(q A v) , q A v) = s((B(q),q)(B(v),v) — 
(B(q),v) 2 ) = eg(q)g(v). After a convenient rcscaling this gives the first case. 
In case (ii) of Theorem 9.2, according to 9.3, 72.^ (g, v; .,.) = </> A g(cj)\ (q A v), .) , 
where the second factor is only defined up to the addition of a scalar multiple 
of <j), but the result of the wedge product is well defined. Equation (10.3) is 
true in particular for any v E ker<f) and it gives (<f>, q)dh\ q A <f> A g(v, .) = for 
any q E TL. Thus we have dh\ q A <fi A £ = for any £ G V* , which implies, if 
dim V > 3, dh\ q A <j) = for any q E TL. This proves that TL is flat. A rescaling 
gives the formulas as written. QED 

11. Conclusion 

It is time to explain how the various statements we established may be used. 
Suppose we wish to study a system (2.1) on a screen TL. If there is a non-trivial 
algebra of first integrals, we must determine it, and as we explained at the 
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beginning of Section 4, it is reasonable to limit ourselves to the first integrals 
which are polynomial in the velocities. 

Laplace [Lap] gave what is maybe the first example of determination of the 
algebra of the first integrals. The system he considered is the Kepler problem 
in space. He devised a systematic research of polynomial first integrals, and 
developed it until he found all the classical first integrals. They are at most 
quadratic in the velocities. Actually Laplace already knew the first integrals he 
found. They all appear in an earlier work by Lagrange [Lag], including the so- 
called Laplace-Runge-Lenz vector. But the method is interesting and general. 
The first step is to determine an a priori information on the higher degree term 
of the integral (compare for example [Whi], p. 332). 

In Section 6 we put this information in a simple and general form: after ho- 
mogenization, the higher order term is a multilinear form, whose symmetry 
corresponds to a rectangular Young tableau with two rows. In other words it is 
a polynomial in q A v, the projective impulsion. 

Suppose we find a quadratic first integral, i.e. a function G = T — U , where U 
is a function of the position only, T a function of the position and the velocity 
which is homogeneous polynomial of degree 2 in the velocity. Such a first integral 
has a lot to do with a Hamiltonian, as emphasized in the case of a flat TL by 
Lundmark. Translated into the language of projective dynamics, Lundmark's 
remark is simply: if G is not the Hamiltonian, it may become the Hamiltonian 
after a change of screen. 

We proposed a simple test to decide if a screen exists that makes G the Hamil- 
tonian. The condition we found is almost exactly the characterization of Lund- 
mark's cofactor systems. Projective dynamics simplifies the statements, gives a 
geometrical interpretation and brings to these theories the powerful tools from 
Young theory (compare [Ben]). 

The test is as follows. We must decide if L = T + U is a Lagrangian for some 
screen. Following Proposition 7.8, we just look at T. We homogenize T in an 
Ra £ 'P'jfiV), i- c - m a quadrilinear form having the same symmetries as the 
Riemann curvature tensor in Riemannian geometry. This symmetry corresponds 
to a 2 x 2 Young tableau. 

The quadrilinear form Ra does not depend on the screen. Guided by Proposition 
10.2 we ask if Ra has a non-trivial kernel. In this case G is not a Hamiltonian, 
but according to Proposition 7.5 it can still give a preserved prc-symplcctic form 
to the system. This form exists if the reduced pair (Ra,, W°) in Proposition 10.2 
is a compatible pair, which we decide using Proposition 10.3. 

If the kernel of Ra is trivial, we use Proposition 10.3 to see if there is a screen 
such that G is a Hamiltonian, and which is the screen. The screen must be 
contained in a hyperplane or a quadric. Proposition 10.3 provides a justification 
for the definition of a cofactor system. 

The theory continues asking if there is another first integral, if it is a Hamiltonian 
for some screen, if with the first it forms a "cofactor pair" . This would imply bi- 
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hamiltonianity, separability and integrability. Thus for a system with n degrees 
of freedom, two quadratic first integrals may be enough to integrate. All this is 
explained in the thesis of Lundmark. 

Acknowledgements. I wish to thank A.V. Borisov, A. Chcncincr, B. Coll, M. 
Garay, I.S. Mamaev, R. Montgomery, M. Przybylska and J.M. Pozo Soler. 
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